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Abstract 

Starting with a nearest-neighbors tight-binding model, we rigorously investigate the bulk 
zero-field orbital susceptibility of a non-interacting Bloch electrons gas in graphene-like solids 
at fixed temperature and density of particles. In the zero-temperature limit and in the semi- 
conducting situation, we derive a complete expression which holds for an arbitrary number 
of bands with possible degeneracies. In the particular case of a two-bands gapped model, all 
involved quantities are explicitly written down. Besides the formula that we obtain have the 
special feature to be suitable for numerical computations since it only involves the eigenvalues 
and associated eigenfunctions of the Bloch Hamiltonian, together with the derivatives (up to 
the second order) w.r.t. the quasi- momentum of the matrix-elements of the Bloch Hamilto- 
nian. Finally we give a simple application for the two-bands gapped model by considering 
the case of a dispersion law which is linear w.r.t. the quasi-momentum in the gapless limit. 
Through this instance, the origin of the singularity, which expresses as a Dirac delta function 
of the Fermi energy, implied by the McClure's formula in monolayer graphene is discussed. 
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1 Introduction and the main results. 



1.1 Introduction. 

Isolated for the first time in 2004 by Novosclov et al. , graphcne is one of the most intriguing 
sohd because of its exceptional electronic properties which, in great part, originate from its purely 
two-dimensional structure; see [3] for a review of these properties. Regarded to be a semimetal 
since the pioneer works of Wallace in 1947 [45] (he was the first one who investigated the band 
structure of graphcne through a nearest-neighbors tight-binding model) , graphene has the feature 
to exhibit a very high diamagnetic susceptibility. Note that this characteristic is shared by most 
semimetals, and to a lesser degree, by most narrow-gap intrinsic semiconductors. 

From a theoretical viewpoint, the first attempts to account for this large diamagnetic behavior 
go back at least to the 50's with the works of Smoluchowski 04] and Hove [21]. Even though 
these papers only were concerned with graphite (viewed as a piling up of graphene sheets) which 
shows the same peculiarity than graphene (graphite is known as one of the strongest diamagnetic 
solids), they are nevertheless interesting in regards to the method. The starting point in [331[3T] is 
the formula of Peierls j361 Eq. (62c)] obtained in the tight-binding approximation (it amounts to 
the same thing to consider a single-band model) for the zero-field orbital susceptibility of conduc- 
tion electrons at non-zero temperature, together with the Wallace band structure |45| (two bands 
touching at the Brillouin zone corners: the lower one completely filled, the upper completely empty 
at zero temperature). However in 1953, the Peierls formula threw back into doubt by Adams [2]. 
This latter pointed out the existence of other contributions coming from the interband effects 
induced by the magnetic field; but apart from particular cases in which these contributions are 
far from being negligible, no general formulas were derived. Note that the quest to a complete 
formula (including the interband contributions) of the zero-field susceptibility of Bloch eletrons 
in solids generated an intense activity up to the middle of 70's, see e.g. [28] for a review. Note 
also that the first rigorous investigations in this direction came as late as 1990, see [19]. Recently, 
Briet et al. [Hj have rigorously derived a complete formula for the zero-field susceptibility of Bloch 
electrons at fixed temperature and density for both metals and semiconductors/insulators. They 
have also given a rigorous justification to the Landau-Peierls approximation in low temperature 
and density limit. 

In 1956, McClure wrote down in the high temperature regime a formula for the zero-field orbital 
susceptibility of Bloch electrons in a two-dimensional honeycomb lattice (modeling a single-layer 
graphite) taking into account the interband effects, see [IS] Eq. (3.15)]. To achieve that, he first 
calculated the energy levels induced by the magnetic field for the two-dimensional Wallace band 
structure in |45j . This band structure yields a dispersion law which is, as a first approximation, 
linear w.r.t. the quasi-momentum near the band-touching points at the Brillouin zone corners. 
Subsequently McClure used these calculations to work out the Helmholtz free energy from which 
he derived the zero-field orbital susceptibility. He stressed that 'all the conduction electrons dia- 
magnetism of the two-dimensional model is due to band-to-band transitions'. 

Following [26], most papers were mainly concentrated either on 'real' graphite, or on com- 
pounds involving graphite- like structure, see e.g. [27l SU |39l |40]. All the same we mention that 
when dealing with (infinite) single- layer graphite model, Safran et al. [39], and later on [40], recov- 
ered the McClure formula in a regime of small values of the chemical potential (they do not tend 
to show some restrictions on the temperature regime). However the question of the behavior of the 
zero-field orbital susceptibility at zero temperature is not addressed. It should be point out that 
the method used in [27] |43l [39l [40] basically is different from the one implemented in [26] . Indeed 
their starting point is a given formula for the zero-field susceptibility of a 3-D Bloch electrons gas, 
in that case the one derived by Fukuyama |17| . Afterwards they carry out some transformations 
in that formula to take account of the dispersion law induced by the band structure of materials. 

In 1987, Levintovitch at al. [25] were interested in the zero-field orbital susceptibility of a 
two-dimensional Bloch electrons gas in narrow- gap semiconductors at zero temperature and fixed 
density. Their method consisted in expressing a (quite) modified version of the zero-field suscepti- 
bility formula derived by Wannier et al. }46j only in terms of a dispersion law corresponding to a 
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two-bands gapped model (the origin of energy is taken at the middle of the gap, and when the gap 
width is zero, the two-bands model used by MeClure in f26| is recovered). In the zero-temperature 
limit together with the semicondueting situation, they found that the zero-field orbital susceptibil- 
ity solely is determined by the interband contributions, and further, it is in inverse proportion to 
the gap width, see [251 Eqs. (22)-(23)]. Besides in the gapless limit, it has a singularity expressed 
as a Dirac delta-like dependence on the Fermi energy (that is the zero-field orbital susceptibility 
is minus infinity at band-touching point, and vanishes otherwise). They claimed that their result 
is in good agreement with the McClure's formula in the zero-temperature limit. 

After the experimental discovery of Novoselov et al. [33] , many investigations arose to figure 
out the large diamagnetic behavior of monolayer graphene. Their approach substantially is dif- 
ferent from the previous authors: it is based on the observation that the dispersion law (derived 
from the 2-D Wallace band structure in [JS]) in the neighborhood of the band-touching points in 
the Brillouin zone is, due to its linearity, analogous to the energy of massless relativistic particles. 
Since that suggests a massless Dirac fermions-like behavior of quasi-particles around the band- 
touching points, most papers were focused on the diamagnetism of massless Dirac fermions, see 
e.g. [22j [m [29l [23]. We do not give further details about their method since this present paper 
will not deal with massless fermions (we anyway mention that the starting point in [18l [29] is the 
zero-field susceptibility formula derived by Fukuyama in |17j). However we emphasize that for a 
gapless Dirac fermions system, Koshino et al. and Fukuyama obtained that the zero-field orbital 
susceptibility at zero temperature is expressed as a delta function of the Fermi energy (the origin 
of energy is taken at band-touching point), see [221 Eq. (37)] and [HI Eq. (3c)]. Besides they 
stated that this result is in accordance with the McClure's results at finite temperature. In [55], 
Nakamura was interested in gapped Dirac fermions systems. In the semiconducting situation, he 
found that the zero-field orbital susceptibility at zero temperature is in inverse proportion to the 
energy gap, see [5S1 Eq. (7)]. In the zero-energy gap limit, the results of [H] [TH] are recovered. 
For the sake of completeness, we mention that within the framework of 2-D systems with a linear, 
massless Dirac-like dispersion law, the behavior of the zero-field susceptibility at zero-temperature 
(brought up just above) for such systems with and without a band gap, was already pointed out 
by Sharapov et al. in 2004, see [H Eqs. (7.4)- (7.6)] and also [41]. 

In the light of the above mentioned results, the singular behavior of the zero- field orbital suscep- 
tibility at zero temperature of Bloch electrons in graphene seems to be regularized when opening 
a gap in the band structure of graphene. This brings us naturally to be interested in the zero-field 
orbital susceptibility of gapped graphene to shed light on the possible origins of the delta-like sin- 
gularity, other than the purely semimetallic nature of graphene at zero temperature. Furthermore, 
the study of fundamental properties of gapped graphene has recently found a renewal of interest 
due to new technical processes allowing to open a band gap in graphene materials, see e.g. |34[|35] . 

In this present paper, we rigorously investigate the zero-field orbital susceptibility of a non- 
interacting Bloch electrons gas in graphcne-like solids. By a graphene-like solid, we mean a purely 
two-dimensional crystal which is an intrinsic semiconductor (possibly with a narrow gap) at zero 
temperature. In the one-body approximation, our starting point is a nearest-neighbors tight- 
binding model for the Hamiltonian of a Bloch electron in a finite-size graphene-like solid plunged 
into a constant magnetic field applied perpendicular to the solid plane. Our model of tight-binding 
Hamiltonian is a Harper-like operator (see e.g. [32l [TTl [15] ) , and is similar to the one of Brynildsen 
et al. [2 who lately have rigorously investigated the Faraday rotation in graphene-like solids. From 
the thermodynamic limit of the grand-canonical pressure (the question of its existence is there 
fully addressed, see Lemma [TTT] below) . we derive a general formula for the bulk zero-field orbital 
susceptibility at fixed positive temperature and density of electrons, see (jl.lip . This is made 
possible by the use of the so-called gauge invariant magnetic perturbation theory (see e.g. [10 ] [31 ] 
and [31 [121 [131 [T] [8] for further applications), followed by the Bloch-Floquet decomposition. After 
carrying out some convenient transformations needed to perform the zero-temperature limit, then 
in the semiconducting situation, we get from (|l.lip a complete formula for the zero-field orbital 
susceptibility at zero temperature and fixed density which holds for an arbitrary number of bands 
with possible degeneracies, see (|1.15p . It consists merely of ()1.16p when considering the particular 
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case of a two-bands gapped model. The whole of terms involved in p.l6p are explicitly written 
down and are expressed only in terms of eigenvalues and associated eigenfunctions of the Bloch 
Hamiltonian, together with the derivatives (up to the second order) w.r.t. the quasi-momentum 
of its matrix-elements. The terms involved in (|1.15p have the same special feature what makes our 
formulas suitable for numerical computations. Finally, we give an application for the two-bands 
gapped model by treating the case of a dispersion law which is linear w.r.t. the quasi-momentum 
in the gapless limit (this is a simple 'textbook model' for monolayer graphene). The singular 
behavior of the bulk zero-field orbital susceptibility at zero temperature obtained in the gapless 
limit is construed from a mathematical viewpoint. 

1.2 The setting and the main results. 

Consider a two-dimensional finite-size crystal subjected to a uniform external magnetic field 
applied orthogonal to the solid. By crystal, we mean an ideal crystalline ordered solid composed 
of a regular lattice of identical atoms whose nuclei are fixed at their equilibrium positions. Due to 
the magnetic field, the Bloch electrons possess an orbital magnetic moment together with a spin 
magnetic moment. Since we only are interested in orbital diamagnetic effects, we disregard the 
electron spin. Furthermore we neglect the self-interactions: the Bloch electrons gas is supposed 
strongly diluted. Besides the finite-size crystal is at equilibrium with a thermal and particles bath. 

Let us precise our assumptions. Let T := {v : v = X]fe=i ''^k^Lk', fik € Z} be a Bravais lattice 
generated by two independent R^-vectors a/j, k = 1,2. Let i7 := {x : |x| < inf.yg-c\{o} 1^ ^ '^^1} 
be the Wigner-Seitz (W-S) cell of T, and denote by Q the finite set of vectors modeling the basis 
(i.e. the positions of the ions at equilibrium in fJ). Hereafter, we suppose that Card(0) > 2. 
The finite-size ion-core mesh is given by An := {x + u : x G 9, i^ = J^t^i ^k^k', \nk\ < N}, 
N G N*. Note that it basically boils down to the same thing to consider a central region of 
A := Q + Y = {x + V : X £ Q,v G T} which stands for the ion-core mesh of the corresponding 
infinite crystal. The magnetic field is given by the 3-dimensional vector B := {0,0, B), B > 0. 
By seeing a R^-vector x = {xi,X2) as the 3-dimensional vector (a;i,a;2,0), the magnetic potential 
vector which reigns in the crystal is defined by Ba.{x) := ^{0,0, B) A (xi,a;2,0) = -y(— X2,xi,0) 
(we use the symmetric transverse gauge). Throughout this paper, we denote by b :~ eB/c G 
the electron cyclotron frequency where we set mdcctron = 1- Note that below we set ft = 1 too. 

Under the above conditions, let us introduce the one-particle Hamiltonian. To fix notation 
and give some additional assumptions, we firstly start with the situation in which the crystal is 
infinite. In view of our discrete approach, let ^^(A) be the one-particle Hilbert space and denote 
by {i5x}xgA the canonical basis, where i5x(y) = 1 when y = x and otherwise. Recall that the 
kernel of any bounded operators A on /^(A) is given by A{x,y) := {Sx,ASy) (x, y) £ A^, where 
(• , •) denotes the usual scalar product on /^(A). In the one-body approximation, let Hq be the 
Hamiltonian on (A) which determines the dynamics of each Bloch electron in the absence of the 
magnetic field. We denote by Ho{- , • ) : A x A — > C the kernel of Hq. In the whole of this paper, 
we assume that Hq satisfies the following three properties (given through its kernel): 

(PI) Hq is a self-adjoint operator, i.e.: 

V(x,y)eA2, i/o(x,y) = i/o(y,x). 

(P2) Hq is a nearest-neighbour tight-binding operator, that is: 

3C>0 : Ha{x,y)^0 when|x-y|>C. (1.1) 

(P3) Hq commutes with the translations of the T Bravais lattice, i.e.: 

VueT, HQ{x + v,y + v) = Ho{x,y) (x,y)eA^ 

Usually the magnitude of the non-zero hopping matrix-elements for Hq{- , • ) is fixed to some finite 
constant (typically, fixed to 1). In this case, and we suppose that hereafter, the kernel Ho{- , • ) is 
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EAD (^Exponentially Almost Diagonal), i.e. there exist two constants ci,C2 > s.t. 

V(x,y) e A\ |i7o(x,y)| < cie-=^l-''l. (1.2) 

In presence of magnetic fields, the Hamiltonian Hi, on P(A) is defined through its kernel iJb(- , • ) : 
A X A — > C which is obtained from Hq{- , • ) by using the so-called 'Peierls substitution': 

V(x,y)eA2, i/,(x,y) :=e""^(^''')i/o(x,y) beR+,h=l- (1.3) 

here <^(- , ■ ) stands for the usual magnetic phase defined by 0(x, y) := x • a(y) = ^{yiX2 — Xiy2)- 
Since the phase is antisymmetric, then Hi, is a self-adjoint operator by virtue of (PI). Moreover 
(|l.ip still holds true for i7h( - , • ), and due to (P3), iJf, commutes with the magnetic translations: 

-iv e T, i?b(x + y + vy'"*'^^'''^ = e''"^(''^"'i7fc(x, y) (x, y) G A^. (1.4) 

Defined via (|1.3p . Hi, actually is a Harper-like operator; see e.g. [32l[Tl] for its spectral properties. 

Let us get back to our initial problem dealing with the Bloch electrons gas in the finite-size 
crystal. The one-particle Hamiltonian in presence of magnetic fields is defined from Hb by: 

HN,b := XNHbXN b e M+, (1.5) 

where xn denotes the characteristic function of the finite-size ion-core mesh Aat. For any b e M^, 
(II. 5p defines a family of bounded and self-adjoint operators on P{A) for any N G N* . Obviously 
this definition corresponds to choose Dirichlet boundary conditions on the edges of A^r . Hereafter 
we denote by {eb^^}jLi the set of eigenvalues of H^.t counting multiplicities and in increasing 
order. 

Now define the pressure of the Bloch electrons gas in the finite-size crystal under the grand- 
canonical conditions. Recall that in the grand-canonical ensemble, the fixed external parameters 
are (/3, z, Card(A7v)), where: /3 := (fc^T)"^ > is the 'inverse' temperature (fcs stands for the 
Boltzmann constant) and z := e*^'' > is the fugacity (/x G K stands for the chemical potential). 
For any f3 > 0, z > and 6 > 0, the grand-canonical pressure is given by (see e.g. [37]): 

Let us introduce some additional notation. Let i?f,(- ,• ;^) : A x A — ?> C be the kernel of the 
resolvent operator Rb{£,) := (i?b — ^'S ^ p{Hb)- Let F be a simple positively oriented closed 
contour around the spectrum of Hb defined as in (|2.ip . We set f(/3, z; • ) := ln(l + ze"^' ) for any 
/3,z>0. 

Our main technical result is given in the following lemma: 

Lemma 1.1. (i). For any /? > 0, z > and b > 0, the thermodynamic limit of the grand- canonical 
pressure exists. More precisely: 

PiP,z,b) :. Jim P.(/3,z,6) = ^^^^ d^ f(/3, z; ^ E i?.(2., 2.; 0, d-T) 

uniformly in (/3, z) G [/3i, /32] x [zi, Z2], < /3i < < 00 and < zi < Z2 < 00. 

(ii) . For any /? > and b > 0, P{P, • , 6) is a C°° -function on M!j_. 

(iii) . For any /? > and z > 0, P(/3, z, • ) is a -function on M-|_. 

By virtue of Lemma 1 1.1 1 introduce the bulk density of particles defined by: 

dP 

p{P,z,b):=pz—{/3,z,b), /3>0, z>0, 6>0, (1.8) 
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as well as the bulk orbital susceptibility defined as the second derivative of the bulk pressure w.r.t. 
the intensity B of the magnetic field (see e.g. [T1[3J[5J[7] and references therein): 

eVd^P 



X{(3,z,b):= {-] -^il3,z,b), /3>0, z>0, 6>0. (1.9) 



From now on, the density of particles po > becomes our fixed external parameter. Seeing the bulk 
density p as a function of the /x- variable (instead of the z-variable), we denote by /x^"-' (/?, po,b) € R 
the unique solution of the equation p(/3,e^'',6) = po- Note that the inversion of the relation 
between the bulk density and the chemical potential is ensured by the fact that p i~> p(/3, e*^'', 6) 
is a strictly increasing function on (0,00). Under these conditions, the bulk zero- field orbital 
susceptibility at fixed positive 'temperature' /? > and density po > is defined by: 

XiP, po) = X{(3, PO, 0) := e^^""^P'Po,b=°\b = 0). 

Before formulating the main results of this paper, let us introduce some additional notation. 
Let T* be the dual lattice of T generated by the vectors hi, I = 1,2 defined from the a^'s by 
afc • b; = 2n5k.i- Denote by f2* the Brillouin zone of T*. By virtue of (P3), the Floquet theory 
for periodic operators allows us to use the bands structure of the spectrum of Hq, see e.g. [24] . 
Under the Bloch-Floquet unitary transformation, P{h.) = J^-^, dkP{Q) and Hq can be seen as the 

direct integral J^, dki7(k), see e.g. [38l Sect. XIII. 16]. For each k G il* , the Bloch Hamiltonian 
i/(k) lives on and it is defined via its kernel which reads as (see e.g. [6l Eq. (2.4)]): 

V(x,y) e i/(x,y;k) := ^ e-^'^-(^+-^)i/o(x + t^, y) kef!*. 

Note that in the canonical basis of ^^(8), the matrix of ff(k) is nothing but [H{xi,x^;'k)]i^m, 
I < l,m< Card(0). Introduce as well on the operators {dk„H)(k) and (9^-^ ) (k) with 

a, 7 = 1,2 and k e O*, generated via their kernels respectively defined on by: 

(afc„i?)(x,y;k) := -^if(x,y;k) = ^^Y1 yJe-*''-(^+"^^) Ho(x + t^, y), (I-IO) 

{dk,d,^H)i^,y;,k) :- ^^^i/(x,y;k) = 

= - E (^7 + "7 - +Vo.- yJc-^'^-(^+-y)i/o(x + v,y) = ^£l^F(x,y; k). 

Since Ho{- ,■) is BAD, the kernels {dk^H){- , • ; k) and {dk^dk^H){- , • ; k) are well-defined on 0^. 

Getting back to the Bloch Hamiltonian, we denote by {Ej(\i)}'^^^^^^^ the set of its eigenvalues, 
counting multiplicities and in increasing order. Due to this choice of labelling, the Bloch ener- 
gies Ej{-) are T*-periodic and continuous functions but they are not differentiable at (possible) 
crossing-points; for a review see e.g. [301 Sect. III]. Hereafter we will consider both situations: 

(Al) The Ej{-)^s, j = 1, . . . , Card(0) are non-degenerate for Lebesgue-almost all k e SI*. 

(A2) The Ej{-)^s, j = 1, . . . , Card(0) are non-degenerate everywhere on ft* . 

In other words, (Al) means that the set formed by the crossing-points has a zero Lebesgue- 
mcasure. As for (A2), it amounts to the same thing to assume that the Bloch bands are simple; 
the j-th Bloch band is defined by £j :~ Uken* ^jO^)- Under either assumption, the spectrum 
of Hq is absolutely continuous and given as a set of points by (j{Hq) = [J^^^'^^^ £j. The energy 
bands of Hq corresponds to the disjoint union of the fj's. Under the condition (Al), the f/s may 
overlap, be contained in each other, or even coincide. Generally speaking if maxfj < minfj-i-i for 
some j — \,. . . , Card(0) — 1, then we have a non-trivial spectral gap. If maxfj = minfj+i, then 
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wc have a trivial spectral gap. 

Lastly, denote by Epipo) the Fermi energy which is defined by £f{pq) '■= lim^^oo m'"^ Po, 0). 
This limit always exists and actually defines an increasing function of poj see [S] Thm. 1.1]. 

The main results of this paper are collected in the following theorem: 
Theorem 1.2. Let po > be fixed. 

(i) . For any P > 0, let p!^^^ = p^^\l3,po,0) G R &e the unique solution of the equation p{f3,e^'^ ,0) = 
Pq . Then the bulk zero- field orbital susceptibility at fixed positive temperature and density can he 
written as: 

where for each k e fl* , denoting by i?(k;^) := (-ff(k) — ^)^^ the fiber of the resolvent: 

V,{k;0 R{k;0{{dk,H){k)R{k;0{dk,H){k) - {dk,H){k)R{k;0{dk,H){k)}x 

X R{k;0{{dk,H){k)R{k;0{dk,H){k) - {dk,H){k)R{k;Oidk,H){k)}R{k;0, (1-12) 

p4(k;C) := R{k;0{{dlH){k)R{k;0{dk,H){k)R{k;0{dk,H){k) + 

- {dk, dk, H) (k)i?(k; [ {dk, H) (k)i?(k; (,) {dk, H) (k) + (9,., H) (k)i?(k; (dk, H) (k)] + 

+ {dl H) (k)i?(k; {dk, H) (k)i?(k; ^ H) (k) } i?(k; ^ , (1.13) 

P3(k;0 := -i?(k;e){^a\i/)(k)i?(k;0(a,2,i?)(k) + \{dlH){k)R{k-i){dlH){k)+ 

- {dk,dk,H){k)R{k-0{dk,dk,H){k)]R{k-i). (1.14) 

(ii) . Suppose (Al). Assume the semiconducting situation, i.e. the Fermi energy lies in the middle 
of a non-trivial gap, that is there exists M g {1, . . . , Card(0) — 1} s.t. max^M < minfAZ+i o,nd 
£f{po) ~ (maxi£M+inin£A/+i)/2. Then there exist 2M functions ) and Oj,o(' ), j — I, . . . , M 
defined on fl* outside a set of Lebesgue-measure zero s.t. 

Xipo):^ lim X{P,po) 

e\'l 1 /• „ ... . , .... (1-15) 



f dkJ2 {5,4 (k) + {E,{k) - £F{po)}dM\, 



cj 4Card(e)|17*| ^ 
and the above integrand can be extended by continuity to the whole of Q* . 

(iii). Suppose (A2) together with Card(0) = 2. Assume the semiconducting situation, i.e. the 
Fermi energy lies in the middle of the spectral gap separating the band £1 from £2 . Then in the 
zero-temperature limit, the bulk zero-field orbital susceptibility at fixed density is given by: 

x( \ (''\^^ 1 I ^.,^r ui/k) , {i;i(k)-gf^(po)}Pi,,(k) ] 

^(^")-Uj 8^L'''g\(i^.(k)-i^.(k))'+^+ (i..(k)-i.,(k))'- j' (^-1^) 

where the functions Ui^;(- ) and ) are respectively defined in p.29p - p.3ip and (j3.32[) - p.34p . 

Remark 1. All functions 5j,i(-) and c)j^o(-), j = 1,...,M appearing in (|1.15p can be explic- 
itly written down. They only involve the eigenvalues and associated eigenfunctions of the Bloch 
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Hamiltonian -ff(k), together with the derivatives (up to the second order) w.r.t. the fee-variables 
(a = 1, 2) of the Card(0) x Card(6)-niatrix elements of -ff (k), see Section 3.2.1. As for the func- 
tions Ui./(-) and Oi_;(-), I = 0,1,2, they have the same peculiarity, see page [H] for their explicit 
expressions. This special feature makes our formulas suitable for numerical computations. 

Remark 2. We mention that the assumption (Al) in (ii) is artificial in the following sense. One 
actually can prove that (jl.lSP still holds true when the Ej^s are degenerate on a subset of ft* 
with full Lebesgue-measure. However in that case the functions ) and fj.o(' ), j = ^, ■ ■ ■ , ^ 
can not be expressed as it is mentioned in Remark 1: their expressions are more complicated and 
require the use of the orthogonal projection corresponding to each Ej. For further details, see [5]. 

Remark 3. When (Al) is replaced with (A2) in (ii), then (jl.lSp still holds true with the same 
functions ), j — 1, . . . , M and / = 0, 1. But in this instance these functions are smooth in the 
k- variable by virtue of the analytic perturbation theory, see e.g. [351 Sect. XII]. In the particular 
case of Card(0) = 2, one has the following identities: 

Remark 4. As it is common practice in physics, the starting point in this kind of problem is a 
given matrix for the Bloch Hamiltonian H(k). Since the trace in the formula (|l.lip only involves 
the derivatives (up to the second order) w.r.t. the fc^-variables (a = 1, 2) of the Card(8) x Card(0)- 
matrix elements of the Bloch Hamiltonian together with its 'inverse matrix', then a such formula 
allows 'hand'-calculations when dealing notably with 2x2 matrices. More precisely, the calcula- 
tions from p. lip are quite easy when considering dispersion laws which are polynomial in k in the 
gapless limit. Under these conditions, the zero-temperature limit simply follows from the formula 
p.ll[) by the application of the residue theorem without resorting to the formula (|1.15|) or (|1.16p . 
see the following paragraph for some examples. Besides we mention that under the assumption 
(Al), we give another formula for the zero-field orbital susceptibility at fixed positive temperature 
and density which have the same special feature than the one mentioned for the formula (|1.15|) , see 
(|3.1ip in Proposition 13.21 In the particular case of Card(0) = 2, all the terms involved in (|3.1ip 
are explicitly written down, see Lemma 13.81 Furthermore the assumption (A2) makes possible 
the rewriting of (|3.1ip as a sum of two terms: the first one is the so-called Peierls contribution, 
the second one stands for the interband contributions, see (|3.38p in Proposition 13.91 In the zero- 
temperature limit and in the semiconducting situation, only a large number of the terms involved 
in the interband contributions will give rise to p. 161) . 

Remark 5. The regularity properties announced in (ii)-(iii) of Lemma 11.11 are far from being 
optimum. On the one hand, one can prove that z i— >■ P(/3, z, b) can be analytically extended to 
the complex domain C \ (— oo, — e^'^^ '] where i?^^'' := inf o'(iJb), see e.g. [31 SI [7]. On the other 
hand, the use of the gauge invariant magnetic perturbation theory to prove (iii) allows us actually 
to get that b i— >■ P{/3, z, b) is a C°°-function. For further details, see the proof of Proposition 12.31 

1.3 A simple application for the two-bands gapped model - Discussions. 

The purpose of this paragraph is to present and discuss some calculations on the bulk zero-field 
orbital susceptibility of Bloch electrons at zero temperature and fixed density in the framework 
of a two-bands gapped model for which the dispersion law is linear w.r.t. the quasi-momcntum in 
the gapless limit. For monolayer graphene, this model is a case in point as a first approximation. 

We first give a series of calculations from two distinct Bloch Hamiltonians, see (|1.17p and (|1.2ip 
below, which both have the same eigenvalues. Subsequently, we discuss the results. 
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Calculation (1). The starting point is the foUowing Bloch Hamihonian: 

^''^^y-=(kj.k. ''^t') ^>o,^^^*. (1.17) 



Let i?i(k) := — y^P + lkp and E2(k) := — £'i(k) be its two eigenvahies. The matrix-resolvent is: 

Note that the second derivatives w.r.t. the fc^-variable, a = 1,2 of the matrix elements in (|1.17p 
are all zero. Therefore, each quantity in (|1.13p and (|1.14l) is reduced to the 2 x 2-zero matrix. 
After carrying out some calculations involving products of 2 x 2-matrices, one gets for the quantity 
defined in (|1.12p (which we denote by T^g'' (k; ^) to express the dependence on the model in ([LT7])): 

^""'^^-"(^^(kFoWky^U-^fc^ ^-s ) "^^"^ (^-i^) 

Let /3 > and po > be fixed. From (|1.18p . the formula leads to: 

Suppose the semiconducting situation, i.e. the Fermi energy lies in the middle of the spectral gap 
(the width of the gap is 2S) separating the band £i from £2- By performing first the integration 
w.r.t. the ^-variable by the residue theorem, and subsequently the zero-temperature limit (for 
further details in regards to the method, see the proof of Theorem 11.21 (ii) — (iii)), one has: 

X{po) := ^hm A'(/3,po) = -^(^) J^J^^ dk{5i,i(k) + {E,ik) - £F{po))^,,o(k)} , 

1 (52 

where: 5ii(k)=: tH and Oio(k)=0. (1.19) 

' ((52 + |k|2)i (^2 + |k|2)f ^'"^ ^ ^ ' 

Finally, by performing first the integration w.r.t. the k-variable (note that the presence of the (S's 
makes integrable the function 5i^i(- ) near k = 0) and only after the gapless limit, one obtains: 

lim A'(po) = -00. (1.20) 
Calculation (2). The starting point is the following diagonal Bloch Hamiltonian: 

(,) _/£;i(k) 0\ r i?i(k):=-V^^W, ^>okGf2* (121) 

The matrix-resolvent is: 

(£;i(k)-e)(i?2(k)-e) I i?i(k)-ej 



Let us remark that: 



^ ' ' Jp + W V 1/ 



Due to this feature, it is easy to see that the quantity in (|1.12p is reduced to the 2 x 2-zero matrix. 
After working out the trace of the quantities in (|1.13p - (|1.14p . (jl.lip yields: 



y|Rn^ MM 1 1^ /• .uA (-i)^W' ( ^ \ f f(/3,M'"^e) 



iS^ + m\2i7rJ Jr ^(i?,(k)-0 
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Let us mention that the apphcation of the residue theorem for each integral w.r.t. f (inside the 
braces) provides us with a derivative of f (/3, /x*^"' , • ) of order three and two respectively. This 
feature leads, in the zero-temperature limit and in the semiconducting situation, to the following: 

^'(po):-!™ <¥(/?, po)=0. (1.22) 

p— foo 

Let us now discuss these results. Starting with the Bloch Hamiltonian in (|1.17|) . the singularity 
expressed as a delta function of the Fermi energy is recovered, see (|1.20p and Section 1.1. Through 
the calculations leading to (|1.22l) . this singularity can be accounted for by the following. 

The Bloch Hamiltonian-matrix in (|1.17p is related to the one in (|1.2ip by: 

H^^\k) = P{k)H^'^\k)P~\k) k e n*, (1.23) 
where the change-of-basis matrix together with its inverse read for instance as: 



P(k) := ( <5+VP+W^ j , F"\k) := 



V2^(52 + |k|2 2^<52 + |k|2 



(1.24) 



Denote by P|''^(k;^) the quantity defined in ([112)) but with H'-'^^k). By replacing i?(')(k) with 
the r.h.s. of (fOSi in V^^^ (k; ^ which only survives in the calculation (1), one gets: 

(k;0= (k;e)P-i(k) + l^f'^)(k;e) ken*. (1.25) 

Here w!^'"'^\k] £,) consists of a sum of terms, each of them containing at least a 'derivative w.r.t. 
the fcc-variable of the change-of-basis matrix' (or its inverse). A generical term is for example: 

P(k)i?('') (k; 0P'\^) (dk, P) ik)H^^^ (k)i?('^) (k; (dk, H^''^ ) (k)i?('*) (k; ^ (dk, H^"^ ) (k) x 

X i?('^)(k;e)i/('')(k)(5fc,p-i)(k)P(k)i?('^)(k;0P-'(k) = 1,2, t ^ j. 

As mentioned in calculation (2), vi,'^\k;^) = (in the matrices-sense). In view of (|1.25p and 
(|1.18p . this means that only Wg''''' (k; ^) (through its trace) provides a non-zero contribution to the 
zero-field orbital susceptibility at zero-temperature, see (|1.19p . Therefore, the singular behavior 
appearing in the gapless limit (|1.20p is closely related to the fact that the maps k i— > P(k) and 
k i—> P^^(k) are non-differentiable in k = when 6 = 0, see ()1.24p . 

Finally, the identically zero result obtained in (|1.22p tends to show that our 'linear model' as 
a first approximation for monolayer graphene is not good enough. Some numerical computations 
based on a more 'realistic' dispersion law derived from a 'gapped two-dimensional Wallace'-type 
band structure is needed. 



1.4 The content of the paper. 

Our current paper is organized as follows. Section 2 is dedicated to the proof of Lemma 11.11 
The proof of (i) essentially is based on a simplified version of the geometric perturbation theory 
recently developed in [13]. Since the method for the discrete case is fully detailed in [H Sect. 3.3], 
the proof of (|1.7p just is outlined. The crucial ingredient involved in the proof of (iii) is the gauge 
invariant magnetic perturbation theory applied to the resolvent operator Rb{£,)- It allows us to 
keep a good control over the linear growth induced by the magnetic vector potential when dealing 
with the diagonal part of the kernel Rb{ - ,■ ;£,)■ Section 3 is devoted to the proof of Theorem 11.21 
In Section 3.1 we focus on the bulk zero-field orbital susceptibility at fixed positive temperature 
and density. In particular, we prove Theorem ll.2l (i) from the results of Lemma ll.ll (iii) together 
with the Bloch-Floquet decomposition. In Sections 3.2 and 3.3 we prove Theorem II. 21 (ii) and (iii) 
respectively. It essentially follows the outline of the proof of [SJ Thm. 1.2 (i)]. From the formula 
(jl.lip . we separately consider both situations corresponding to each one of the assumptions (Al)- 
(A2) and we perform the zero-temperature limit in the semiconducting situation. 
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2 Proof of Lemma 11.11 



2.1 Proof of (i) - (ii). 

Under the grand-eanonieal conditions, let /3 :— {ksT)^^ > and z := e*^^ > 0. 
For any & > 0, let F be the eounter-clockwise oriented simple closed contour defined by: 

T:^me[S^,S+],-s^^±^}um = S±, cj^e[_^,^]}, (2.1) 

where (5_ := inf cr(_ffh) — 1, :~ supcr(_fffc) + 1. By construction F surrounds the spectrum of 
Hi,, and the closed subset surrounding by F is a strict subset of S) := G C : 3C G (— 7r//3, tt//?)} 
which is the holomorphic domain of the function ^ i— >■ f(/3, z; <^) :— ln(l + ze^'^^) for any z > 0. 
We point out the fact that the reals S± in p.ip can actually be chosen 6-independent since the 
Z^-norm of H), is bounded from above by some constant uniformly in b. This follows from the Schur- 
Holmgren criterion in Eq. (1.1)] applied to the kernel (|1.3I) . combined with the estimate (|1.2p . 
Note also that under our conditions, there exists a constant c = c(/3, z) > s.t. 

VeeF, |f(/3,z;OI<ce-''^«. (2.2) 

Now we use the Dunford functional calculus |16[ Sect. VI. 3] to write down a more convenient 
formula for the grand-canonical pressure of the Bloch electrons gas in the finite-size crystal. 

First of all we need to give some properties on the resolvent operator of -ff^v.b defined in 
P31) . Let iV G N* and 5 > 0. For any ^ € p{HM,b), denote by RN,b{.C) — {Hn^b - ■ 
Obviously RN,b{0 lives on the Hilbert space l^(Ajv), and moreover it owns a kernel denoted 
by RN,b{','',0 ■ X Ajv —J- C. This kernel is EAD (=Exponentially Almost Diagonal) in 
the following sense: for any compact subset K of p{HN,b), there exist two constants ci,C2 > 
independent of b,N s.t. 

V(x,y)eA^, sup|i?^,,(x,y;e)| <cie-^^l^-''L (2.3) 

The estimate (|2.3p simply follows from the fact that the kernel of H^^b is EAD, see [H Prop. 3.5]. 
An important ingredient for the following is the uniform estimate on the trace of i?7v,&(0- 

sup|Trp(A„)(i?w.6(0)| <c, (2.4) 



for some if-dependent constant c > 0. To obtain ()2.4p . all we have to do is write the trace as the 
sum over A^r of the diagonal matrix-elements of i?jv^b(^), and afterwards use the estimate (j2.3|) . 
Let /3 > 0, z > 0, 6 > and N G N*. With F the contour as in ((2?T|) . introduce on P{An): 



The Dunford functional calculus provides us with the identification £Ar(/3, z, b) = ln(lAr-|-ze '3^" '') 
which holds on P{Ais[). In view of (|1.6p . this allows us to define the grand-canonical pressure by: 

P.(/3,.,^)4^^^Tr(£.(/3,z,6)).i^^^ /dCf(/3,z;0Tr(i?.,(C)). (2.5) 

Clearly the trace in the r.h.s. of the first equality makes sense (do not forget that we in fact deal 
with N X A^-matrices). The commutation of the complex integral w.r.t. ^ with the trace in the 
r.h.s. of the second equality is ensured by the Tonelli's theorem, see estimates (|2.4p and (|2.2p . 
We draw the attention to the fact that if the following limit holds: 

lim \ . Tr(i?jv,b(e))- ^ !,^n^ E^fefe^'^)' (^.6) 
jv->oo Lard(AAr) Card(fc)) 
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then (|1.7I) (only in the pointwise-sense) wiU follow from the dominated convergence theorem. 

Let us outline the proof of (|2.6p which uses a simplified version of the geometric perturbation 
theory in [T31 Sect. 4]. All the details about this method in the discrete case can be found in [51 
Sect. 3.3.1]. Below [-J denotes the floor function. Pick e S (0,1), and assume that N is large 

enough. The first step consists in dividing An into a 'core part' An and an 'edge part' An, where: 

2 

Aat :={x + t^ :xe e,t^ = ^nfeafc;|nfe| < iV - [A^^J}, In Ajv \ A^f. 

fc=i 

From a geometrical argument, one can see that An contains (2(A^ — [A^*^]) + 1)^ W-S cells. This 
implies that the number of W-S cells in A^v is of order 0{N^^'^) for A^ sufficiently large. Owing to 
this splitting of An, we expect only the contribution to Ti{RN,b{0) coming from the core region 
A^v to give rise to the limit in (j2.6p . With this aim in view, the geometric perturbation method 
consists in approximating the resolvent RN,biO with the operator UN.biO defined by: 

UnAO Xa.MOXa, + Xx^RnAOxx^^ (2.7) 

where and are the characteristic functions of Ajv and Ajv respectively. Afterwards, it 
remains to control the behavior when A^ oo of the trace of the r.h.s. of (|2.7p . On the one hand: 

^^^^Tr(x = ^i?^,.(OXAj = (2iV+l)'card(e) ^ ^^^.^(-'-5^) = ^(]^)' 

xSAiv 

where we used (|2.3p together with the fact that the number of W-S cells in Aat is of order 0(A^^+'). 
On the other hand, by using that Rb{£,) commutes with the magnetic translations, see (jl.4[) : 

Card(A^) (2Ar-M)2Card(e) ^^^"^-'"'^^ " C^rf(e) "^-'"'^^ 

To complete the proof of (|2.6p . we need to control the behavior of the trace of {RnAO ~ UN,b{0}- 
One can prove that (Card(AAr))"iTr(i?Ar,b(^)-J77v,&(0) = 0{N-^). We do not give further details. 

In order to complete the proof of Lemma [1.11 (i). consider for any (/?, z) € [/3i,/32] x [21,22], 
with < /3i < /32 < 00 and < < Z2 < 00, the following quantity: 

Q(/3,2,6):=-^ / def(/3,z;g)f J. ^ i?^,,(x, x; ~ „ \,^. Y.^'>{^^^-^o), 
271 J Y VCard(Ajv) Card(e) ^ J 

where the F-contour is defined as in ()2.ip but with j32 instead of (3. Then V(/3, z) e [/3i , /32] x [zi , Z2] : 

\QiP,z,b)\<c [ |de|e-^^^« — -i-— i?^,b(x,x;e)- „ ^ ^"fe^' ^) 

Jr Card(Ajv) ^f^^ Card(e) ^ 

for some constant c = c(/3i, 22) > 0. By virtue of (|2.6p . the proof of Lemma [TTT] (i) now is over. 

Let us conclude this paragraph by proving (ii) of Lemma 11.11 Let /? > be fixed. For any 
^ G 2) (the holomorphic domain of f), z h-> f(/3, z; ^) is a C°°-function on (0, 00). In particular: 

and ^ I— > (9™f)(/3, z; 5), m G N* is holomorphic on T). Moreover it obeys the following estimate: 

VCeF, |(an)(/3,^;e)l <ce-"''^« meW, 

for some constant c = c(m) > 0. The proof of Lemma 1 1.1 1 (ii) follows by standard arguments. 
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2.2 Proof of (iii). 

The proof of (iii) essentially is based on the application of the so-called gauge invariant magnetic 
perturbation theory to the resolvent operator Rb{£,)j see Lemma [2.21 below. This method makes 
it possible not only to investigate the regularity of the kernel i?f,(- , • ; ^) w.r.t. the 6- variable, but 
also to write down formulas for its partial derivatives w.r.t. &, see Proposition 12.31 For reader's 
convenience we collect in the appendix of this section all the proofs of intermediary results needed 
to prove (iii). 

For the sake of simplicity, let us see 6 as a parameter on the whole real line. Pick a 6o G K. 
For any ^ G p{Hbo) and & G M, introduce on P(A) the operators i?f,(^) and Tf,(^) through their 
kernel which respectively are defined by: 

V(x,y)GA2, 4(x,y;e) :-c'*^*(^''')i?bo(x,y;e) 6b:^b-bo, (2.8) 
f,(x,y;0 e'*'"^^-''') ^{^.^MCx.z.y) _ i}Hb„i^,z)Rb„{z,y;0, (2-9) 

zGA 

where, for any arbitrary vectors u, v, w G A, fl(u, v, w) :— (j){u, v) + (^(v, w) + </'(w, u). Note that 
fl(u, V, w) stands for the flux of a unit magnetic field through the triangle generated by u, v, w. 
From (|2.8p . RbiO clearly is bounded on P(A) since its kernel is BAD. Ditto for Tb{^): 

Lemma 2.1. Let fep G K fee fixed. For each compact subset K C p{Hi,^^) there exists > s.t.: 

V6gM, sup ||fb(e)|| < c^l^fel. (2.10) 

The main point of the gauge invariant magnetic perturbation theory is the below identity 
(|2.1ip . It asserts that for any b sufficiently close to bo kept fixed, the resolvent Rb{0 <^an be 
approximated on /^(A) by Rb{(,)', the /^-norm of the 'corrective term' being of order ©(jfe — feo|): 

Lemma 2.2. Let bo be fixed. Then for each compact subset K C p{Hb„) there exists > 
s.t. in the bounded operators sense on /^(A); 

VfeG [feo-?x,6o + ^K], V^G 7^, RbiO = RbiO - RbiOTbiO- (2.11) 

The above identity, twice iterated and written in the kernels sense, allows us to establish: 

Proposition 2.3. Let 6o G M and K C p{Hbg) be a compact subset. Then there exists > s.t. 
b I— >■ i?;,(x, y; ^) is a -function on {bo — <j/<-, feo + '^k) G K and V(x, y) G A^. In particular, its 
second partial derivative at bo reads as: 

V(x, y) G A^ {d^,R)b„ (x, y; = - ^^^^Y^i?^„ (x, y; + 

+ X! ^feo(xiZi;^){(^(x,zi) + (/i(zi,y)}fl(zi,Z2,y)-ff6o(zi,Z2)i?bo(z2,y;C) + 

zi,Z2£A 

X] ■^''0 ^1 ' ' =22,23)^^60 (zi , 22)^60 (22,23; 0^(23, Z4, y)Hbo (23, Z4)i?6„ (Z4, y; £_)+ 

Zl ,...,Z4gA 

+ ^ H i?b„(x,zi;e){fl(zi,Z2,y)}2i7f,„(zi,Z2)i?6„(z2,y;0, (2-12) 

Zl ,Z2GA 

and the first both terms in the r.h.s. of p.l2|) vanish when x = y fey antisymmetry of the phase. 

Let us prove Lemma [1.11 (iii). Let /3 > 0, z > and feg G M. With the contour F defined in 
(|2.ip (with bo in place of 6), recall that the bulk pressure of the Bloch electrons gas reads as: 
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From Proposition [531 'wc know that there exists a <;^r > such that T encloses the spectrum of Hi, 
for all b satisfying |6 — 5o| < "Tr, and moreover b i-^ i?b(x, x; ^) is a C^-function on (60 — Tr, f'o + "Tr) 
e r and Vx e A. On the other hand, from the expressions (|2.12p and (|2.15p (with y = x) and 
by using the same arguments than the ones used in the proof of Lemma |2 . II below . then one has: 

Vxe A, V6e (5o-<rr,6o + <^r), sup |(9^i?)b(x, x; <e)| < c / = 0,1,2, 

for some F-dependent constant c > 0. Lemma 1 1.11 (iii) now follows by standard arguments. 

2.3 Appendix: proofs of intermediate results. 
Proof of Lemma 12.11 From the definition (|2.9p . it holds on A^: 

|ff,(x, y; 01 < \Sb\ ^ |x - z||iffc„ (x, z)| |z - y\\Rt„ (z, y; 01, 

where we used that |c"' — 1| < Vx e M, followed by the estimate |fl(x, z,y)| < |x — z||z — y|. 
Now use that i?bo(' > ' ) is EAD, and therefore Rtgi' r \ too by [5J Prop. 3.5], to get ride of the 
factor |x — z| and |z — y| respectively. Ergo there exist two X-depcndent constants ci, C2 > s.t.: 

Vee A', \niK,y;0\<c,\Sb\c-'-\^-y\. 
The lemma follows by the Shur-Holmgren criterion in [iTl Eq. (1.1)]. ■ 
Proof of Lemma l2.2l Let K C p{Hi,g) be a compact subset. Pick a £^ E K. From (|2.8p . one has: 
V(x, y) e A^, {HbMm^, y) = E e'*^(^('''^)+^(^'''»e^''°^('''^)i/o(x, z)i?,„ (z, y; 0- 

Note that 0(x, z) + 0(z, y) = fi(x, z, y) + (f>{x, y). Then from (fOj) together with (|2^ : 

(Hfci?b(0)(x, y) - e'^^^(--y) ^ i/fc„(x, z)i?b, (z, y; + fb(x, y; e). 

zeA 

Since (i7f,oi?t„ (0)(x, y) = (5x,y + ?-Rho(x,y) with (5x,y := ^xi,yi<^3;2,y2 (the Kronecker symbol), then: 

V(x, y) e A^, {{H, - OMO)i^, y) = e^""^(-'y)5x,y + r6(x, y; 0- 

This implies in the bounded operators sense on Z^(A): 

{Hb-OMO^i + TbiO- (2.13) 

Next use that there exists a constant > s.t. K C p{Hb) for all \b — 6o| £ Cx, see fTTJ Thm 1.1 
(i)]. Thus for such 6's, we can invert the operator (Hi,—^) with G K as above. The proof is over.B 

Proof of Proposition [2T3l Let 5o £ M and K C p{Hbo) be a compact subset. Pick a ^ £ K. 
The starting point is the identity (|2.1ip twice iterated. Written in the kernels sense, one has 
V6 e M satisfying | & - 5ol < '^a' and V(x, y) G A^: 

i?b(x,y;0 = i?b(x,y;e) - (i?6(C)Tfc(0)(x, y) + [MOTKm^.Y) - (^^^(OTfc (0)(x, y). (2.14) 

For each kernel in the r.h.s. of p.l4p . we expand the exponential phase factor appearing in the 
kernels of RbiS.) and Tf,(^) in Taylor series up to the second order in Sb. Then we have: 

i?, (x, y ; = E {Sb)'^-^^^^^Rb„ (x, y ; ^ + 7^i°) (x, y ; ^ , 

fc=0 
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{MOniOK^, y) = E (x, zi; e){«<56 - [Sbf [(0(x, zi) + 0(zi , y)) + -fl(zi, Z2, y)] } x 
zi,z2eA 

X fl(zi,Z2,y)F6o(zi,Z2)i?6„(z2,y;0 +'^i^'(x,y;^), 
(i?bm2(0)(x,y) = -(<5&)2 5^ i?b„(x,Zi;C)fl(zi,Z2,Z3)iyfc„(zi,Z2)i?6„(z2,Z3;0x 

Zi , . . .,Z4 £A 

X fl(z3,Z4,y)i76„(z3,Z4)i?6o(z4,y;^) +7^[^'(x,y;^), 

where the remainder terms 7?,['^ (x, y; / = 0, 1, 2 satisfy the property that their first two deriva- 
tives at bo are identically zero. Note that we can disregard the last term in the r.h.s. of (|2.14p 
since by virtue of (|2.10p . its ^^-norm behaves like Od^fop) when |i56| — >■ uniformly in ^ € i^. 
Next from the above expansions, for b G [bo — <;k, bo + <iK] sufficiently close to bo, it holds on A^: 

^6(x, y; - Rbo (x, y; = ^^{^^(x, y)Rba (x, y; ^)+ 

-i ^ i?6o(x,zi;^)fl(zi,Z2,y)iJho(zi,Z2)-R6o(z2,y;0} + o('5^)- 

Zl ,Z2 £ A 

Perform the limit b ^ bo, and we get that the map b i— > i?h(x,y;^) is differentiable at bo with: 
{dbR)bo{^, y;0 := *'/'(x,y)i?6o(x,y;0 - * X! ^bo(x,zi;Ofi(zi,z2,y)-ff6o(zi,z2)i?6o(z2,y;0- 

Zi.zoGA 

(2.15) 

This result can be extended on the whole of (^o^'^^a'i ^o+'^^a')- The proposition follows by induction. 
Note that from (|2.1ip iterated n-times and written in the kernels sense, our method allows to get 
the expansion of the kernel of the resolvent in power series of Sb up to the rt-th order. By an 
iterating procedure, one therefore can prove that b i?b(x, y; ^) is a C°°-function near &o, G K 
and V(x, y) e A^. Although the formula obtained for the n-th derivative (9^i?)h„(x, y; ^) is quite 
complicated, its diagonal part only involves the flux fl with the kernels Ho{- , • ) and Ro{- ,■ ; £,). ■ 

3 Proof of Theorem 11.2 . 
3.1 Proof of (i). 

We start by writing down a general formula for the bulk orbital susceptibility under the grand- 
canonical conditions. From (|1.9p together with the results of Section 2.3, one directly obtains: 

XW,z,b)= (-) 1 \ ^ I def(/3,z;g)^(agit;)b(x,x;g) /3 > 0, z > 0, & > 0, (3.1) 
\cj /3 Card(e) 27r 7r ^ 

where {d^R)b{- ,■ ',0 equals (|2.12l) (only the two last terms survive for the diagonal part). 

From now on, the density of particles po > becomes our fixed external parameter. Under the 
conditions of Lemma II. 1[ the relation between the fugacity and the bulk density in (jl.Sp can be 
inverted. This is made possible since for any /? > and b > 0, the map z i— >■ p(/3, z, b) is strictly 
increasing on (0,oo), and actually defines a C°°-diffeomorphism of (0,oo) into itself. Then there 
exists a unique z^^\l3, po,b) € (0, oo), and therefore a unique fj,^°\(3 , po , b) G R satisfying: 

p{(3,e^'^'°'^^'P^'''\b)=po>0. 
In this way, the bulk orbital susceptibility at fixed positive temperature and density is defined by: 

X{P, PO, b) := X{(3, {(3, po,b),b)= X{(3, e^^^'"' ^^'""'"^ ,b) (3>0,b>0. (3.2) 
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Under the zero magnetic field assumption, let us now prepare the formula p.2[) (knowing p.ip ) 
for the Bloch-Floquet decomposition. To do that, introduce the 'position operator' in P{A) as: 

Vx = (a;i,a;2) € A, (X • ei^)(x) = .t,^(x) yj e l^A), I = 1,2, (3.3) 

where lc{A.) denotes the space of /^-functions with compact support. 

For any ^ G p{Hq) and any q, r,s,t £ {1, 2}, define on P(A) the following families of operators: 

Sg.r.sAO ■■= i?o(0[X • e„ [X • e,,i/o]][X • e„ [X • e^, i?o(0]], (3-4) 
Tg.r^sAO ■■= ^o(0[X • e„ Ho][X ■ er, i?o(e)][X • e„i/o][X • e*, RoiO]- (3.5) 

Here [• , • ] denotes the commutator on /^(A). Since the kernels of Hq and i?o(0 ^"^^ EAD, then 
all the above commutators are well-defined on ^^(A). As we will see hereafter, the families of op- 
erators defined in p.4p and (|3.5p have the feature to commute with the translations of the Y-lattice. 

Involving p.4p and (|3.5p . the bulk zero- field orbital susceptibility at fixed positive temperature 
and density can be written as: 

Lemma 3.1. Let /3 > and po > be fixed. Let /i'^^ = /i(°^(/3, ^ = 0) G M 6e the unique 
solution of the equation p{(3,e^'^,b = 0) = Pq. Then one has: 

- 5l,2,l,2(0 - r2,l,2,l(C) + + Tl,2,2,l(0 - ri,2,l,2(0}(x, X)- (3-6) 

Proof. By inserting (|2.12p into p.ip . then by setting 6 = we get: 

X X! I 5Z ■Ro(x,Zi;0{fl(zl,Z2,x)}^^fo(zi,Z2)i?o(z2,X;0 + 

- 2 ^ i?o(x, zi;^)fi(zi,Z2,Z3)i7o(zi,Z2)i?o(z2,Z3;C)fl(z3,Z4,x)i/o(z3,Z4)i?o(z4,x;0 [■ 
zi,...,z4eA ^ 

From the above expression, it remains to give both ingredients which lead to (j3.6p . The first one 
consists in rewriting the magnetic fiux fi(x,y, z), for arbitrary x, y,z G A'^, as the following: 

fl(x,y,z) := 0(x,y) + (f>{y,z) +(p{z,x) = ^[{yi - Zi){x2 - y-z) - {V2 - Z2){xi - yi)], 

where we used that 0(u, v) = \{u2Vi — U1V2) = \[{u2 — V2)vi — {ui — Vi)v2]- The second one 
consists in using the following identity. Let A{- , • ) be the kernel of a bounded operator A. Then: 

V(x,y)GA2, (xi-y,)^(x,y) = [X.e;,A](x,y) / = 1,2, 

where X-ej is the position operator defined in p.3p . Its kernel reads on A^ as (X-e;)(x, y) = aJiJx.y, 
/ = 1,2. The remainder of the proof only is rearranging of terms. Note that 2,1, 2(0 '^^ p.6p 
can be replaced with 5i,2,2,i(0 oi" '^2, 14, 2(0 oi" '^2,1,2,1(0 due to the above arguments. ■ 

It turns out that (|3.6p actually is suitable for the Bloch-Floquet decomposition which is at the 
root of the formula (|l.lip . Let us recall the framework of it, see e.g. [38l Sect. XIII. 16]. 
Let n* be the Brillouin zone of the dual lattice T*. Introduce the 'Bloch isometry' as: 

U : ll{K) ^ ^ := r dkl^ie) 

Jn* 

(Z^V)(x;k) = -^^e-*'^-(^+-)V'(x + t.) kGl7*, xGe, iJElliA), 
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which can be extended by continuity in a unitary operator on P{A). The adjoint of U reads as: 
{U*(t>){^ + v) = -^ j dke''"(2+"V(x;k) xe0, -i^eT, 0(-;k)e/2(0)^ 

Since i/o commutes with the translations of T, then the unitary transformation of Hq is dccom- 
posable into a direct integral UHqU* ~ j^, dki7(k). For each k S fl* , the fiber Hamiltonian 
H{k.) acts on P{Q) and it is defined via its kernel which reads as, see e.g. [51 Eq. (2.4)]: 

V(x,y) e e^, if(x,y;k) := ^ c-'"^-(^+--x)i/o(x + y) k e Q* . 

Note that for each k e il*, the fiber i?(k) is a self-adjoint operator on ^^(6) due to the identity 
ffo(x,y) = Ho{y,x.) (see assumption (PI)) which implies the one ff(x, y;k) = iJ(y,x;k). Here- 
after we denote by i?(k;^) := {H{k) - ^)"^ £ p{Ho) the fiber of the resolvent on P{Q). 

We now can start the actual proof of Theorem 11.21 (i). In view of (|3.6p . we need an expression 
for the fiber of the operator V<j,r,s.t(C)i V = 5 or T. Below we use the shorthand notation Qq 
to denote Hq or i?o(C): C ^ p{Ho)- Due to the assumption (P3), then [X • e^jQc] with a = 1,2, 
also commute with the lattice translations. By induction, ditto for [X • e^,[X • ea,(5o]] with 
q;,7 = 1,2. Therefore the unitary transformation of Vq_r,s.t{C} is decomposable into a direct 
miegraXUVq^r^sACll^* = In* dk Vg,^,^,t(k; C). For each k € fl* the fiber Vg,^,,,t(k; 0, V = 5 or T 
lives on and it is defined by: 

Sg.r^s.t(k;0 ■■= ^(k;0([X • e„ [X • e,, Ho]]){k){[X ■ e„ [X • e*, i?o(0]])(k), 

Tg^r,sA^; ^(k; e)([X • e„ Ho]) (k) ([X • e„ R^iO]) (k)([X • e,. Ho]) (k)([X • e^RoiO]) (k). 

Here the fibers ([X-e^, Qo]){k) and ([X-e-y, [X-e^, Qo]]) (k) with Qo = ^^o or i?o(C) and a,^ = 1, 2 
are defined through their kernel which read for each k G il* and any (x, y) £ respectively as: 



( [X • e„ , go] ) (x, y ; k) : = ^ e-'^'^- (2+"^^) ( [X • e„ , Qo] ) (x + , y ) , 
([X • e„ [X • e„, Qo]])(x, y; k) := ^ e-*''-(^+"-^) ([X • e^, [X • e„, Qo]])(x + v,y). 



Let Qo( - , ■ ) be the kernel of Qo- Since ([X • e^, (5o])(x + v,y) = {x^ + Va — 2/_^)Qo(x + y) and 
([X • e^, [X • eQ,(5o]])(x + '*^,y) = (i^^ + - y^)(x^ + - y^)Qo(x + i^, y), then one has on 6^: 

([X • e„,go])(x,y;k) = i — Q(x,y;k), ([X • e^, [X • e„, Qo]])(x, y; k) = (i— )i — 0(x, y; k), 
where: V(x,y) £ 9^, Q(x,y;k) := ^ e-''' (2+"-y)Oo(x + y) k e 

Here these identifications are made possible since the kernels of Ho and RoiO ^^'^ EAD. Denoting 
by (9fe^Q)(k) and (9fe^9fc^g)(k) the operators on ^^(6) generated respectively by the kernels: 

(9A^.Q)(x,y;k) := — Q(x,y;k), (5fe,9fc„Q)(x,y;k) := ^^-^g(x,y;k) - ^^-^ Q (x, y ; k) , 

then for each k e il* the fiber V<j.r,s.t(k; ^), V = 5 or T can be rewritten as: 

S,^r,sA^\£.) = R{k-0{dk,dkMk){dk^dk,R){k-0, (3.7) 
r,,...,t(k;0 = i?(k;0(9fc,iJ)(k)(5fc,i?)(k;0(afe,iJ)(k)(afc,i?)(k;e). (3.8) 
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Since the kernel of Vg.r,s.t(0 is related to the one of the fiber Vq.r,s.t(k; ^) by: 

V(x,y) e A^, (V,,.,.,t(e))(x,y) = dkc^'^-(''-'')(V,,.,.,t(k; 0)(x, y) V = 5 or T, 

where we have used |14[ Eq. (4.1)], then the expression in p.6p can be rewritten as: 

^(^'^-«)40'^c^^d(W^i 

+ ^5i,i,2,2(k; -51,2,1,2(1^; - r2,i,2,i(k; +r2,i,i,2(k; +Ti,2,2,i(k; ~ 'ri,2,i,2(k; e)}(x, x)- 
Next it remains to insert in p.7p and p.8p these identities: 

(afc„i?)(k;e) = -i?(k;e)(afc„if)(k)i?(k;e) ke 1]*, a = 1,2, 
(9fe^afc„i?)(k;0 = i?(k;0(afc,i?)(k)i?(k;0(5fe„i?)(k)i?(k;0 - R{k;Oidk,dk^H){k)R{k;0+ 
+ R{k; (afc„ H) (k)i?(k; (9fc, ) (k)i?(k; ^) k € 17* , a, 7 = 1 , 2. 

Thus (|1.11[) follows by gathering together the terms by decreasing number of fiber of resolvent. 

3.2 Proof of (ii) - (iii). 

The starting point is the expression of the bulk zero-field orbital susceptibility at fixed positive 
temperature and density obtained in (jl.lip . In view of (jl.lSp we need to subject this formula 
to some transformations again so as to perform the integration w.r.t. the ^-variable. Only after 
we will perform the zero-temperature limit. Let us remark that the presence of fiber resolvents 
in (|l.lip indicates the possibility to bring into play, through their kernel, the eigenvalues and 
associated eigenfunctions of the fiber Hamiltonian. As we will see below, this will turn out to be 
essential to perform the integral w.r.t. the ^-variable by the residue theorem. 

Before beginning, let us fix notation. For each k G fi*, denote by {£'j(k)}^^'^'^^^'' the set of 
eigenvalues of the fiber Hamiltonian i?(k) counting multiplicities and in increasing order. Due 
to this choice of indexation, the Ej^s are T*-periodic and continuous functions, but they are not 
differentiable at (possible) crossing-points. Denote by {uj{- ; k)}^^'^"^'^^ the set of associated eigen- 
functions. They form a complete ortlionormal system in Unlike the Ej^s, the eigenfunctions 
k I— > Uj{- ; k) may be not continuous at (possible) crossing-points since they arc defined up to a 
k-dependent phase factor which cannot be always chosen to be continuous at crossing points. 

For any integer l,m = 1, . . . , Card(0) and a, 7 = 1, 2, let us define the quantities: 

^i,™(«;k) := "/(x;k)[(afe„i?)(k)«„(x;k)] kef!*, (3.9) 
xee 

CT/,„,(a,7;k) := ^^lKxn^Pfc.9fe^if)(k)M™(x;k)] ke^!*. (3.10) 

Due to the phase factors presence in the u/s, both above quantities may be not continuous at 
(possible) crossing points. Besides from p.9p . one has ■ki,m{a\'k) = 7r„j,;(Q;;k) since: 



XI X! QT-^fey;k)w,„(y;k) = X ■ti™(y;k)X — i/(y,x; k)wi(x; k), 



where we used that if (x, y; k) = i7(y, x; k) V(x, y) G Q^- The same property holds for ct;.™ (a, 7; k). 

The sequel of this paragraph is divided into two parts. In the first one, we treat the generical 
situation in which the Ej{- )'s are non-degenerate outside a subset of fi* with Lebesgue- measure 
zero. In the second one, we restrict the number of sites in G to two sites and we consider the 
situation in which both Ei{-) and i?2(- ) are non-degenerate everywhere on fi*. 
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3.2.1 Proof of (ii). 

In what follows we tacitly suppose the assumption (Al), see pageO 
Here is the last rewriting of the zero-field orbital susceptibility before the zero-temperature limit: 

Proposition 3.2. Let /? > and po > be fixed. Let ^'"^ = p^^^/S , po , 0) S R fee the unique 
solution of the equation p{(3,e^^,0) — po. Then for each integer j = 1, . . . , Card(0) there exist 
four families of functions c)j,;(- ) with Z = 0, 1,2,3 defined on fi* outside a set of Lebesgue-measure 
zero, s.t. the integrand below can be extended by continuity to the whole of Q* , and: 

1 / \ 2 , , Card(e) 3 

A-(,,,„,0).-j(f) jj-p^- g y^ dkg^(ftM™;B,(k)fe(k), (3,11) 

with the convention {dlf){l3, p^°'> ; Ej{k)) f{l3, p^°'> ; Ej(k)) = ln(l + e^('^""--^^('^»). 

The special feature of this formula lies in the fact that each function ) can be only ex- 
pressed in terms of eigenvalues and associated eigenfunctions of the Bloch Hamiltonian i?(k), 
together with the derivatives (up to the second order) w.r.t. the fc^-variables {a = 1,2) of 
the Card(0) x Card(0)-matrix elements of H(k). By way of example, the functions i)j.3(-), 
j = 1, . . . , Card(6) are identified in (|3.22p . Note that all functions I = 0,1,2 can also 

be written down, but their explicit expression is not necessary for the proof of Theorem 11.21 (ii) . 
Further, the expansion obtained in (|3.11l) makes the choice of phase for the Uj{- ; k)'s irrelevant. 

Now we give the proof of ProDOsition l3.21 it follows the outline of the proof of [5l Thm 3.1]. For 
reader's convenience we collect in the appendix of this section the proofs of intermediary results. 
For the sake of simplicity, let us rewrite the zero-field orbital susceptibility formula in (|l.lip as: 

2 



where: WK/S, (J<°>, 0) := ^ / (If f(/!,M™;f) / dkTrp|e,(P,(k;{)) ( = 3,4,5. 

The first step consists in expressing the kernel of each fiber resolvent appearing in 7'/(k;^), 
I — 3,4,5 (see formulas (|1.14p - (|1.12p ) in terms of eigenvalues and eigenfunctions of the fiber 
Hamiltonian (via the spectral theorem). Thus the quantities yVi{P,p^^\0) can be rewritten as: 

Lemma 3.3. Let /? > and po > be fixed. Then one has: 

Caid(e) 

W5(/3,/^(°\0) = - / dkC,,,,,,,3,,,(k)x 



ii,... j4=i 



Card(e) 

W4(/3,/^(°\0) = - Yl / dkC,„,,,,3(k)x 



^J_f.. f(/3:M^°^0 ,3.4^ 

2znJrUEM-0'iEM-0{E,A^)~0' ^' ' 



Card(e) 

miM = - ^ E ^ X.akc..„.(k)-/d, (3.15) 
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where Vk G fi*, the functions Cj-^j^Ja.j-ii' ); Cj1.j2.j3i' ) '^'^'^ Cji,j2i' ) '^'"^ respectively defined by: 

''ji J2,i3,i4 (k) := {'^ji j2 (1; k)7rj2 j3 (2; k) — tt^j ^3 (2; k)7rj2 ^3 (1; k)}x 

X {^J3j4(2;k)7r,,j,(l;k) - 7r,3 (1; k)7r,, (2; k)}, (3.16) 

Cji j2 J3 (k) := o-ji.j, (1, 1; k)7rj, ^3 (2; Vi)Ttj^,j, (2; k) + aj.j^ (2, 2; k)7rj, ^3(1; k)7rj3 (1; k) + 

~ ^ju]2{^i 2; k)7rj, j3(l; k)^, ji(2; k) - (Tji ^2(1, 2; k)Trj2 ^3(2; k)7Tj3 ^-,(1; k), (3.17) 

C,,,,,(k) := -5R{a,,.j,(l,l;k)a,,,,,(2,2;k)} + a,,j,(l,2;k)<Tj,.j,(l,2;^ (3.18) 

The quantities 7ri^m(a; k) and ai^mic^, 7; k), a, 7 = 1,2 are respectively defined in p.9p OTid (|3.10p . 

The second step consists in applying the residue theorem in p.l3p - p.l5p . It straightforwardly 
provides us with: 

Lemma 3.4. For each integer ji = 1, . . . , Card(0), there exist four families of functions ) 
with Z = 0, 1, 2, 3 defined on il* outside a set of Lebesgue-measure zero s.t.: 

Card(e) 

W5(/3,A^^°\0) = - V / dk V^(/3,/.W.£;^.^(k))a^.^^,(k). (3.19) 



E / dk^-^(/3,/.(°);£;,,(k))a,,,(k). 



/n particular, for Lebesgue-almost all k G fJ* the function aji,3(- ) is defined by: 

Card(e) 



J_ \ ^ ^ji ji j2 ji (k) 



J2 = 
j2#il 



Lemma 3.5. For each integer ji = 1, . . . , Card(0), there exist four families of functions &ji,i(- ) 
with Z = 0, 1, 2, 3 defined on £7* outside a set of Lebesgue-measure zero s.t.: 

Card(e) 3 

W4(/3,A^(°),0) = - 5] / dk^-i(/3,/.(°);i?,,(k))b,,,(k). (3.20) 
/n particular, for (Lebesgue-almost) all k G fJ* the function bji,3(- ) is defined by: 



f'ji,3(k) := ^Cjijiji(k). 



Lemma 3.6. For each integer ji = 1, . . . , Card(0), there exist three families of functions Cj-^^i{- ) 
with 1 — 0,1,2 defined on f2* outside a set of Lebesgue-measure zero s.t.: 

Card(e) 2 

W3(A/i("\0) = - J2 / dk^-L(/3,^(o);i?,.,(k))c,,.Kk). (3.21) 

We stress the point that all functions appearing in the expansions p.l9[) - (|3.2ip can be written 
down since a such result only is based on identities provided by the residue theorem. This will be 
done in the second part of this paragraph when considering the case of Card(0) = 2. 

Thus Lemmas 13. 4[ and [XBl provide us with an expansion of the type announced in p.lip . 
where the coefficients read for ji = 1, . . . , Card(9) and almost everywhere on f2* as: 

,Card(e) , 



32 
j2=^j 



0,,,,(k) :=a,,,,(k) + b,,,/(k) + c,,,z(k) ? = 0,1,2. (3.23) 
To conclude the proof of Proposition 13.21 we need the following result: 
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Lemma 3.7. The functions ciji./(- ), bji,;(- ) with Z = 0, 1, 2, 3 and Cjj.;(- ), / = 0, 1, 2 are bounded 
and continuous on any compact subset of f2* where Ej-^ is isolated from the rest of the spectrum. 

We do not give a proof of Lemma 13.71 since all details can be found in the proof of [H Lem. 
3.7]. Thereby the functions ) in (|3.22p - (|3.23p might be singular on a set of Lebesgue-measure 

zero where Ej-^ can touch the neighboring bands. However the whole integrand in p. lip comes 
from the complex integrals p.l3p - p.l5p which have no local singularities in k (they actually are 
bounded uniformly in k, see the proof of Lemma l3.3p . The proof of Proposition 13.21 is now over. 

We now can start the actual proof of Theorem 11.21 (ii); it follows the outline of Sect. 4.1]. 
Let po > be fixed. Let = n^°^l3,po,0) e R be the unique solution of p{l3,el^'' ,0) = po- Let 

fFDW,f^^°^;0 ■= (e'^^^"''*'"^ + l)"^ = -/3"H5cf)(/3:M^"^;0 be the Fermi-Dirac (F-D) distribution. 
By involving the F-D distribution and its derivatives, the expansion p. lip can be rewritten as: 

-1 / \ 2 , Card(e) ^ 2 ^ 

-^(^'^-O^^iiy C^^dM^ g y^/k|g^(/3,,(0);i.,(k)).,,,,(kH 

-if(/3,p(°);i?,(k))0,,o(k)|. (3.24) 

Consider the semiconducting situation, i.e. assume that the Fermi energy lies in the middle of 
a non-trivial gap. Then there exists Me {1, . . . , Card(6) — 1} such that max^M < min^M-i-i, 
and moreover, lim^_>oo /Lt^''-'(/3, po) = (max^M +min£M+i)/2. From the pointwise convergences 
O Eqs. (4.2)-(4.3)] together with the estimates (SJ Eq. (4.4)], then all terms in p.24p containing 
derivatives of the Fermi-Dirac distribution will converge to zero in the limit /3 — > oo. This leads 
to (frT5|) . 



3.2.2 Proof of (iii). 

In this paragraph, wc implicitly suppose the assumption (A2), see page El 
Supposing (A2) instead of (Al) does not bring anything to the expansion obtained in p. lip , 
with the only difference that the functions j = 1, . . . , Card(9) and / = 0,1,2,3 (which 

are the same) are in this instance smooth w.r.t. the k-variable. That is ensured by the analytic 
perturbation theory, see e.g. [38l Sect. XII] and also [30l Thm 3.5]. Nevertheless, considering 
the particular case of Card(0) = 2 allows us to write down concise and explicit formulas for the 
functions appearing in an expansion of type p. lip . Here is the counterpart of Proposition 13.21 
when Card(e) = 2: 

Lemma 3.8. Let /? > and po > be fixed. Let p'-^'> = p^^\(3,pQ,0) £ R be the unique solution 
of the equation p(/3,e'^^,0) = po. Then: 

'^(APo,0) = -iQ'^i^|^^dk^g(/3,p(°);i?,(k))6,,(k^ (3.25) 
where Vk G fl* , the functions Oj,;(- ), j = 1, 2, Z = 0, 1, 2, 3 are given by (below m = 1, 2, to 7^ j): 



" /i\ _ 1 f Cj,j^m,m(k) + Cj^m,j,m(k) + Cj^TO,rre,j(k) ^ Cj.m, j (k) + Cj j-.m (k) ^ ^ \ 

21 iEM-E,ik)y + E^{k)-E,ik) 

(3.27) 

^-^("'^ = S(i?„(k)%(k))"+i' - £(£;„,(k)%(k))«+- ^'-''^ 
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and the functions Uj^„(- ) and tij_„(- ), j ~ 1,2, n ~ 0, 1, 2 are given by (below m = 1, 2, m 7^ j): 

u,-o(k) :==C,,™(k), (3.29) 
Uj.i(k) :=Cj^TO.m(k) +Cmjj(k) -Cj>„j(k) -Cjj>„(k), (3.30) 

^j,2(k) ^^j,m,m,m (k) C^y^j- j j(k) Cj ^^(k) ^^j;,?n,j,7n(k) ^_;,m,Tn j' (k), (3.31) 



Oj-o(k)~0, (3.32) 
Oj,i(k) := 2Cj ,ji,m(k) + 2Cm,j,j(k) Cj j(k) C^j- ,fj(k) Cm,j,m(k) ,ji,j(k), (3.33) 

Oj, 2(k) 2Cj,„,™,™(k) - 2C™jj- j(k). (3.34) 

T/ie functions Cj-^j^JaJi (' ); ^ji j2 (' ) '^^^ (• ), ji, 32, jz, ji = 1,2 are defined in Lemma UTSl 

The above lemma is obtained from Lemma 13.31 by setting Card(0) = 2, then by mimicking 
the proof of Lemmas 13.4113.61 Since this is simply calculations based on identities provided by 
the residue theorem, we do not give further details. Note that the above functions with 
j = 1, 2, / = 0, 1, 2, 3 exactly have the same peculiarity than the one mentioned for the functions 
5j,;(- ) appearing in the expansion p.lip . see Proposition 13.21 

The non-degeneracy assumption for the Ej^s, j = 1, 2 allows us to use the regular perturbation 
theory, which from the Hellman-Feynman formula (see e.g. [20l Thm 4.1]), provides the identities: 



^jj(a;k) a =1,2, (3.35) 



- "^■^-("' c.; k) + 2 )^ _ a = 1,2, (3.36) 

d'E,{\,) ^ jR{7r,„,,-(a;k)7r,-,.„(7;k)} 9^g,-(k) ^„ 

^^ = a,,(a,7,k) + 2X, e,{\.) - E^ ^ «,7-l,2. (3.37) 

The derivatives of higher order can also be identified, but we will not use them in the following. 

The use of the regular perturbation theory plays a crucial role in the next result: it leads to 
the identification of the interband contributions coming from the expansion p.25p : 

Proposition 3.9. Let /? > and po > be fixed. Let /i'"^ = p^^^/S , po , 0) e M fee the unique 
solution of the equation p{/3,c^^ ,0) = po. Then p.25p can be split into two contributions: 

X{I3, Po, 0) = A'p(/3, Po, 0) + A'ib(/?, Po, 0), (3.38) 
with, by using the shorthand notations ^(/3, p'*'^; • ) = I = 0,1,2; 

I feV 1 1^ f d^f ^^.J d^Ej{k)d'^E^{k) fd'^Ejik) 



dkl dkl \ dkidk 



A'ib(/3,po,0):=-^(-j _-^^^^dk|g(£;,(k))3,(k)+^g(£;,(k))5,,(k)^, (3.39) 



1 / e \ 1 1 

i 

where Vk e U,* , the functions 3j(- ), j = 1,2 are given by (below to = 1, 2 with m ^ j) 



^^^''^ ■= 2^1 imr^M? + 'E,ik)-E,Ak)j^ ^^-^'^ 



and the functions £'j,;(- ), with j = 1,2 and / = 0, 1 are defined in Lemma \3.8\ 
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The proof of the above resuU can be found in the appendix of this section. The term Xp{(3, po, 0) 
is nothing but the so-called Peierls contribution, see e.g. [311 Eq. (62c)] and Eq. (4.10)]. As 
for the term Aib(/3, po, 0), it stands for the interband (or band-to-band) contributions. As we will 
see below, only the functions multiplying {d^^f){f3, fi^^^ ; ■ ) with I = 0,1 (and therefore, only a great 
part of these interband contributions) will give rise to (|1.16l) . 

We now can start the actual proof of Theorem ll.2l (iii). By using the F-D distribution together 
with its partial derivatives, (|3.25|) reads as: 

Xid, PO, 0) = ^ (^) E I, { E M^"'; (k))5,,+i (k) + 

j — 1 /— 

-^f(/3,/.(°);i?,(k))3,,o(k)|. (3.41) 

Assume that the Fermi energy is in the middle of the gap separating the band £i from £2 ■ From [SJ 
Eqs. (4.2)-(4.3)] together with [51 Eq. (4.4)] again, then all terms in (|3.4ip containing derivatives 
of the Fermi-Dirac distribution will converge to zero in the limit /3 — > 00. Hence, we get (|1.16p . 
Note that the expansion in p.38p also leads to (|1.16p in the semiconducting situation. 



3.3 Appendix: proofs of the intermediate results. 

Proof of Lemma 13.31 We start with the quantity W5 (/3, /x^^-* , 0) . Take a generical term of it: 

W5{f3,p^'\0) := ^^dCf(/3,/^(°);0^^/kTrp(e)(p5(k;0), with: 

P5(k; ■■= ^(k; 0(5fc„i/)(k)i?(k; 0(5fc,i/)(k)i?(k; 0{dk.H){k)R{k; 0{dk^H){k)R{k; C), 

where 0,7 = 1,2. Without loss of generality, let us take for the moment 7 = (needed to simplify 
notations). By the mean of kernels, the trace of p,5(k; ^) can be written as: 

E E n E (a;c„i?)(z2„_i,Z2,„;k)i?(z2„„Z2,„+i;k,0)i?(z8,x;k,C). 

(3.42) 

Next use that the kernel of i?(k;^) can be expressed in terms of eigenfunetions of H{k) as: 



Card(e) 



„ x-^ ■", (x: k)u, (y: k) 

V(x,y)ee^ i?(x,y;k,0= E %{k) ^ ^ ^ ^ ^* ' 

The replacement of each kernel of the resolvent in p.42p with the above expression yields: 



=<5, 



Card(e) 



4 

n E "i"^fe2m-i;k) E (<9fc„-f^)(z2m-l>Z2m;k)Uj„ + i(z2„j;k). 



4 

X 

771 — 1 Z2^^_^ Z^,„e0 



From the definition p.9p . one therefore obtains: 

Card(e) . , 

(3.43) 
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where Cjj^ j2.j3.jj (k) := 7rjj.j2(a; k)7rj2.j3 (7; k)7rj3.j.j (7; k)7rjj.jj (a; k), k G fl* . Let us now prove: 

f(/3,/^(°);e) 



il>i2J3J4(l^) 



< C, (3.44) 



for some k-independent constant C > 0, what ahows to change the order of integration in p.43p . 
We only give the two main ingredients leading to (|3.44p . Let ^0 > so that —^0 < inf o'{Ho) aiid 
large enough. On the one hand, from [T31 Eq. (4.8)] there exists a C = C^o(/3,/x^°^) > s.t.: 



lf(/3,M^°^;e)l 

\E,Ak)-e\E,,ik)-C\---\E,,{k)-^\ 

On the other hand, one has the following rough estimate. There exists another C = C^g > s.t.: 

\h,m{a; k)| < C{Ei {k) + Co) /, m = 1, . . . , Card(e), a = 1, 2. 

This follows from [HJ Eq. (4.6)] together with the fact that (9fe^i?)(k) is bounded from 
to ?°°(0). To see that, use that \\{dk^H){k)\\2^^ = esssup^ge X^yee \i(^kaH){x,y;k)\'^ combined 
with (|1.10p and the estimate (|1.2p . It remains to involve the continuity of the Ej{- )'s to get p.44p . 
To conclude the proof for the quantity W^ifH, 0), it is sufficient to use that it is made up of a 
linear combination of ?i;5(/3, /x^*'' , 0)-like terms. The quantities WiifH, ij.^^\0), with ^ = 3,4 can be 
treated with similar arguments. ■ 

Proof of Lemma 13.41 Denote the integrand appearing in p.l3p by: 

(0) .^ f(/3,A^'°^;C) . . . . ^ 

At first sight, flji.ja j3 j4 (/3, m'"'; ' ) can have poles from the first order up to at most fifth order (in 
the case when ji = j2 — jz — ^'4). Therefore we expect the integral w.r.t. ^ of 0j1.j2.j3 j4 (/?: m'"'; ' ) 
to make appear partial derivatives of f (/3, /i^''-' ; • ) with order at most equal to four. However the 
factor multiplying (9|f)(/3, /i'*'^; • ) is identically zero. Indeed, in view of the function Cjjj2.j3.jj(- ) 
defined in p.l6p . it is identically zero for the combinations of subscripts: 

ji = J2 = is = ii-, ji = k = J3 7^ J4; ji = J3 = ii 7^ .72 ■ (3.45) 

Ergo the expansion of (|3.13p consists of partial derivatives of f(/3,/^^°); • ) of order at most equal 
to three. Then by virtue of p.45p . the quadruple summation in the r.h.s. of (j3.13p is reduced to: 

Card(e) Card(e) 



f(/3 /i'^^^'O 

i3 5^ji 

Card(e) 



W5(/3,/.(°\0) = - EE/ dkC,„,„,3,,i(k)-i- /dC- 



{E,Ak) - £,nE,,{k) - 



E / dkC,,,,2j3..4(k);^ / dC 
, Jn* 2i7r Jr 

Jl, . . . ,J4 = 1 



1 f(/?,/^(°);C) 



P " {E,, (k) - 0^ (£;,2 (k) - (i?.3 (k) - C) (Ej, (k) - ■ 



at most 3 equal 
subscripts 

By applying the residue theorem in the first term of the above r.h.s., one has: 



Card{e) ^ r i(R (0) n Card(e) 

i3#ji j3#ii 



r(/3, Aft"); i;ji(k))+ others terms involving — y(/3,/x(°V) with Z < 2 



.3!i?,3(k)-i?,,(k)9C^^"'" ' " 9e 
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The function Cj-^j-^j^j^{-) appearing in front of (5|f)(/3, /x^*'^; iS^^ (k)) corresponds to 3!ajj_3(- 
The remainder of the proof only is a plain computation using the residue theorem. 



Proof of Lemma 13.51 Denote the integrand appearing in p.l4p by: 



Ji,j2,j3 e N*. 



Note that tij1.j2.j3 iPi m'^"*! ' ) can have poles from the first order up to at most fourth order (in the 
case when ji = j2 = J3). Therefore we expect the integral w.r.t. ^ of f^jj^.j^ jg (/3, /i^^' ; • ) to make 
appear partial derivatives of f(/3, fJ,^^^; ■ ) with order at most equal to three. 
The triple summation in the r.h.s. of p.l4p can be rewritten: 



Card(e) 

W4(/3,/^(°),0) = - 

ii=i 

Card(e) 

- E 

^ V ' 

at most 2 equal 
subscripts 



dkCji.jiji(k) 



2iTT 



2zn {E,, (k) - O'iE,, (k) - OiEj, (k) - C) ' 



The application of the residue theorem in the first term of the above r.h.s. yields: 



f(/3,A^(0);C) _ 1 d^f 



2t7TJr ^{E,,{k)-0^ 3!9e 
This is only the one term which provides a third-order partial derivative of f(/3, • ). 

Proof of Lemma 13. 6i Denote the integrand appearing in p.isp by: 



Note that ij^.j^ • ) can have poles from the first order up to at most third order (in the 
case when ji ~ ^2). Therefore the integral w.r.t. f of ijijad^jM^^'*; ' ) '^dl make appear partial 
derivatives of f (/3, /i^"-* ; • ) with order at most equal to two. We do not give further details. ■ 

Proof of Proposition [3T9I Consider the quantity 6j.3(k), j ~ 1,2 defined in p.26p . Let m —1,2 
with m 7^ j. From p.l6p and (|3.17l) together with the identities (|3.35p - (|3.37p we have: 



E„,{k)-E,ik) 

'dEj{k)Yl 
2 



dEj(k )Yl 
dki 



d'^EjCk) 



dk2 



(7„(l,l;k) 



a,,,(2,2;k)- 
d'^Ej (k) ] ( dEj {\^)\( dEj (k) 



dkl 



Cjjj(k) =(Tjj(l,l;k) 



gg.(k) 



dki 



%,(2,2;k) 



dkl 

gg.(k) 
dkx 



dk2 



'^j,.;-(l:2;k) 



d^Ej{\^) 



-2a,,,(l,2;k) 



dkidk2 
dE,{k)\fdE,{k) 



dkl 



dk. 



Thus the quantity 5j^3(k) can be rewritten 0j.3(k) = mj^3(k) +nj.3(k), k G fi* where respectively: 

^ nA ^-Uf dEd^Wd'E.jk) / dE,{k) y d^E,{k) J dE,ik) \f dE,{k) \ d^E,ik) 
n^j,3W- ^2\^^ dkl J dkl \ dk2 ) dkl \ dkl )\ dk2 ) dkidk2 



nj,3(k) := 



1\ { dEj{k) 
dkl 



'^w(2,2;k) 



dE,{k) 
dko 
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Hereafter we use the shorthand notations f(- ) = f(/3, /x^"^; • ). Let us now prove the foUowing: 



dkl dkl 

From [5] Eq. (4.15)] and after rearranging of terms, one has: 



dkidk2 



(3.46) 



£dkg(£;,(k))m,,3(k) = \ |^^dk0(£;,(k)) 



dEj (k) \ / dEj (k) \ d'^Ej (k) 



dkl I \ dk2 I dkidk2 



d^-Ej{\<i)d^Ej{\<i) OEjik.) d d^Ej(k) dEj(k) d d^Ej(k) 



dkl 



dkl dk2 dkidk2 



dk2 dkl dk2dki 



The above identity corresponds to |5J Eq. (4.16)] added with [SJ Eq. (4.13)] by setting N = j. By 
fohowing the same method as the one concluding the proof of [SJ Prop. 4.2], one obtains p.46p . 
Next, let us consider p.27p . In view of p.40p what we still have to do is to prove that: 



(k)) 



E^{k)-E,{k) 



,{k) - E„,{k)' 

On the one hand, from the definition p.l7p together with the identities p.35p - p.37p : 



(3.47) 



EM-E,{k) 



o-„(l,l;k) 



2 dkl 



a,-,(2,2;k) 



h n ^■k^ i^!^^) 



<7,-,(l,2;k) 



'^w(l:2;k) 



d^Ejjk) 
dkidk2 



Then from the definition p.lSp we get after rearranging of terms: 

Cj",m,j(k) +Cjj-_m(k) 



1 f a^f 



E,n{k)~Ej{k) 



C,.,(k)} = i|^/kg(i?,(k)); 



dkf 



-'^.j(2,2;k) 



d'E,{k) 



dkl "^'^ 



a,-,(l,l;k)-2 



d^Ej jk) 
dkidk2 



a,,,(l,2;k)-2 



Ernik)-Ej{k) 



(3.48) 



On the other hand, the definition of nj,3(k) together with the use of [5l Eq. (4.15)] leads to: 

/^^dk|i(i.,(k))n,,3(k).4£dk|i(i.,(k))/ ' f'^^^^y 



d (dE,{k) 



dk2 V dk2 



'^w(ia;k) - 



d (dEj{k) 



dki\ dkl ^"(2'^^^) 

d (dE,{k) 



dkl \ dk2 



'^«(l,2;k) - 



dk2 V dkl 



'^«(l,2;k) 



(3.49) 



By adding p.48p with p.49p and using the definition of Cj,,n,j(k), the l.h.s. of p.47p is equal to: 



dk0(£;,(k)) 



9i;j(k) /a,-„(2,2;k)7r™,,(l;k) 1 d 



dkl 



Ej{k)-EM 



<3-i,m(l,2;k)^„j(2;k) 1 d ^ n o 



E^{k)~E„,ik) 



2,,/..(2,2;k)+ 

dE,{k) ( (l,l;k)7r,„,,(2;k) 
dk2 \ E^{k)-E„,ik) 



Id (l,2;k)7r„,,j(l;k) 19. o i 

-— a„(l, l;k) - -—^^^——^ + __,^.^.(l,2;k)^ 
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To get (|3.47|) . it remains to use this identity (derived from the regular perturbation theory): 
with T,-j(a,a,7;k) -.^ {uj{- ik), {dk^dk^dk^H){k.)uj{- ik)) = T,-,j(a, 7, a; k) = rj-j(7, a, a; k). 
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